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We study the v = 1/2 Chern-Simons system and consider a self-consistent field theory of the 
Singwi-Sjolander type which goes beyond the random phase approximation (RPA). By considering 
the Heisenberg equation of motion for the longitudinal momentum operator, we are able to show 
that the zero-frequency density-density response function vanishes linearly in long wavelength limit 
independent of any approximation. From this analysis, we derive a consistency condition for a 
decoupling of the equal time density-density and density-momentum correlation functions. By 
using the Heisenberg equation of motion of the Wigner distribution function with a decoupling 
of the correlation functions which respects this consistency condition, we calculate the response 
functions of the v = 1/2 system. In our scheme, we get a density-density response function which 
vanishes linearly in the Coulomb case for zero-frequency in the long wavelength limit. Furthermore, 
we derive the compressibility, and the Landau energy as well as the Coulomb energy. These energies 
are in better agreement with numerical and exact results, respectively, than the energies calculated 
in the RPA. 

PACS numbers: 71.10.Pm, 73.43.-f, 71.27.+a 



I. INTRODUCTION 

The combination of an electronic interaction and a 
strong magnetic field in a two-dimensional electron sys- 
tem yields a rich variety of phases. These are best clas- 
sified by the filling factor v, which is the electron density 
divided by the density of a completely filled Landau level. 
In the case of v = 1/2, the behavior of the system resem- 
bles that of a Fermi liquid in the absence of a magnetic 
field, or at small magnetic fields. This effect can be ex- 
plained with a new sort of quasi-particles: at v = 1/2, 
each electron combines with two flux quanta of the mag- 
netic field to form a composite fermion; these composite 
fermions then move in an effective magnetic field which 
is zero on the average. The interpretation of many ex- 
periments supports this picture. We mention transport 
experiments with anti- dots, in which features of the 
resistivity are related to closed loops of the composite 
fermions around the dots Q, and also focusing experi- 
ments |2| . An overview over further experiments can be 
found in (D . A field theoretical formulation of this com- 
posite fermion picture was first established by Halperin, 
Lee, and Read (HLR) [Q as well as Kalmeyer and Zhang 
H . They formulated the Hamiltonian in terms of Chern- 
Simons (CS) transformed electrons and studied within 
the random-phase approximation (RPA) many physical 
quantities. Besides the theory of HLR there are other al- 
ternative formulations of the composite fermionic picture 
which are mainly based on a gauge transformation of the 
CS Hamiltonian f§. 

We used in the theory of HLR for a determina- 
tion of the ground state energy of the v = 1/2 system 
in RPA. There we found an infrared diverging Landau 
energy. This problem was solved in [§ |, Q by taking 



into account the correct normal ordering of the opera- 
tors. We obtained a Landau as well as a Coulomb energy 
in RPA which is in satisfactory agreement with the ex- 
act and numerical results, respectively ||, [To). In Jll| , 
Conti and Chakraborty (CC) improved the calculation 
of the Coulomb energy by taking into account dynam- 
ical correlations through the formalism of Singwi, Tosi, 
Land and Sjolander |l^] (STLS), known as the Singwi- 
Sjolander theory and established first in the calculation of 
the structure factor of the Coulomb theory. This method 
is a generalization of the RPA. In comparison to the RPA 
this method results in a Coulomb energy which is a better 
approximation to the Coulomb energy obtained earlier by 
numerical simulations of interacting electrons in the low- 
est Landau level by Morf and d'Ambrumenil 1 13 and by 
Girlich |Q. CC did not calculate the Landau energy of 
the v — 1/2 system which would also be a very inter- 
esting quantity to be compared with the RPA as well as 
the true result to c /2, the Landau energy per electron. In 
their theory the resulting zero-frequency density-density 
response function vanishes as the square of the wave vec- 
tor in the long wavelength limit. This is in contradic- 
tion to the RPA result where it vanishes linearly for the 
Coulomb interaction In their paper, CC mentioned 
that the behaviour of the zero-frequency response func- 
tion of their theory is similar to the zero-frequency re- 
sponse functions of the alternative formulations of the 
CS theory ^ cited above. But later on, Halperin et al. 
|]l5f showed that the zero-frequency response function of 
these theories also vanishes linearly in the wave vector 
for the Coulomb interaction if one expands the approxi- 
mation to the RPA. 

In this paper, we will show that the quadratic be- 
haviour of the zero-frequency response function in the 
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theory of CC results from the decoupling of the equal 
time density-density and density-momentum correlation 
function. In CS-theories, this needs a careful treatment 
as will be shown by considering the Heisenberg equa- 
tion of motion for the longitudinal momentum opera- 
tor. With the help of this equation, we will get a zero- 
frequency response function of the v = 1/2 system which 
vanishes linearly for the Coulomb interaction in the long- 
wavelength limit. It can be also learned from this proce- 
dure which relation in the decoupling approximations of 
density-density and density-momentum correlation func- 
tions in CS theories should be fulfilled. In the STLS the- 
ory it is the Heisenberg equation of motion of the Wigner 
distribution function which is used to get the response 
functions of the system. In CS-theories, one has to de- 
couple simultaneously a density-density correlation func- 
tion as well as a density-momentum correlation function. 
We will show that the consistency relation for the decou- 
pling of these two functions are not fulfilled in the theory 
of CC. This is the reason for the quadratic behaviour 
of the long wavelength zero-frequency density-density re- 
sponse function in their theory. We will suggest a STLS 
type theory of the CS system, in which the decouplings 
respect the consistency condition. With the help of this 
decoupling, we will calculate the density-density response 
function, the compressibility, the dynamic structure fac- 
tor, the static structure factor and the Coulomb energy, 
as well as the Landau energy. 

The paper is organized as follows: We will introduce in 
section II the CS-Hamiltonian and calculate the Heisen- 
berg equation of motion for the longitudinal momentum. 
Then we get from this equation the static density-density 
response function in the long wavelength limit, and dis- 
cuss the consistency relation for a decoupling of the equal 
time density-density and density-momentum correlation 
functions. We will study in section III the equation of 
motion of the Wigner distribution function and suggest 
a decoupling which respects this relation. 



II. STATIC RESPONSE FUNCTION OF THE 
v = 1/2 CS SYSTEM 

At first, we review the main steps of the HLR approach 
of the v — 1/2 system. Here, we follow in the main the 
notation of CC pl| . The CS transformation for spinless 
fermions is defined by B 



(r) exp 



2i / dr 'arg(r — r ') p(r ') 



(1) 



where ^t(f) is the electron creation operator, ^ + (r) is 
the transformed fermions, p(f) is the density operator of 
the fermion operators , and arg(r) is the angle that r 
forms with the x-axis. The kinetic part of the Hamilto- 
nian is given after the transformation as 



where rrib is the electron band mass and 



8Ai(f) = / dr'Ur- ?) (p{r r ) - p ) . 



(3) 



is the fluctuation of the CS vector potential, po is 
the mean density of the v = 1/2 system and <f>{r) — 
2 Varg(r) = 2 e z x r/r 2 . We used the convention h = 1 
and c = 1 in the above formula (J2j) . By expanding the 
Hamiltonian in (|J) and keeping only terms up to second 
order in the density fluctuations, one gets in the momen- 
tum space 



H 



ata r + 

2m- k k ^ 



- x vr(fc) 1 p(-k): 



(v (k)+v 2 (k)):p(k)p(-k): 



(4) 



Here n{k) is the Fourier transformed momentum oper- 
ator 7?(r) = —i^f(r) + V^(r), at creates a CS-fermion 

with momentum k and vo(k) — 2ire 2 /k is the Coulomb 
interaction where e 2 = e 2 /e. e is the charge of the elec- 
tron and e is the dielectric constant of the background. 
Vi(k) = An/k and v 2 (k) = (47r) 2 p / (rribk 2 ) are CS po- 
tentials. We denote by : A : the normal ordering of the 
operator A. With the help of the CS transformation of 
the electronic Hamiltonian, we thus get a CS Hamilto- 
nian which does not contain a magnetic field. 

First, we will discuss the asymptotics for small wave 
vectors of the zero-frequency density-density response 
function by a method which allows us to study the de- 
coupling of correlation functions in a Singwi-Sjolander 
theory of the half filled lowest Landau level. Singwi and 
Tosi used in [jl6| a method to obtain the compressibility 
sum rule for the Coulomb system from the Heisenberg 
equation of motion of the longitudinal momentum oper- 
ator q/q ■ Tt(q). The commutator CH of q/q ■ Tt(q) with 
the Hamiltonian (0) is 



1 \ ^9 1 

CH(g,i) = V — at , (t)a s ~,Jt) 



TOfc q k+q/2 



(5) 



(M<l')+MQ')) ■ p{q',t)p{q-q',t) 



q q 



rn b \q' J \ q 



ttl(g') 

mb 



^xn(q + q',t)^ p(-q',t): 



+ [- q *{<t + q', *)j p(-q', t) ■■ x q 



-ffkin = 



2m,f, 



2 r*+(r) (-if + 5A(r)j %(r). (2) 



We consider in addition a coupling to an external poten- 
tial V oxt (q, t) . Since we want to study the adiabatic limit 
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we neglect the time derivative in the Heisenberg equation 
of motion 



£ V ext (q, t)p(q- q', t) fii ) = -CR(q, t) . (6) 

q> ^ 



+0(q°){p(q,t)) . 

By inserting this equation in (Q) we get for the static 
linear response in the limit q — > 



We take the expectation value of this equation of motion 
with respect to the ground state of the system (ground 
state for V cxt = 0) and get the following equation which 
is valid to linear order in V cxt 



V ext (q,t)po = -(CR(q,t)), 



1 



(7) 



p (v (q) + v 2 (q)) (p(q,t)) - ip ^- X (n(q,t)) 

nib \q 



(■) c is the cumulant part of the expectation value 
(AB) C = (AB) - (A)(B) where A, B are the operators 7? 
or p. 

In the following we will discuss the leading order in q 
of the terms in (^). By considering only the linear order 
in V ext we get for q — > the following relation (compare 

H) 



(O) 



* (0> V-=o 
opo 



(pM)>, 



(8) 



for an operator O. Here (-)^ is given by the expec- 
tation value for V cxt = 0. We obtain that (CR(q,t)) c /q 
has the asymptotic behaviour O(q )(p(q)) for q — > 0. For 
deriving this asymptotics we take the q — > limit in ev- 
ery cumulant expectation value in (^). Every additive 
term in ([5]) contains either a term proportional to q or 
has a term linear in q' ■ q/q (the terms with a quadratic 
q' in the numerator cancels). We now discuss the lead- 
ing terms in (|t|). We begin with the third term of the 
right hand side of equation (Q) and retransform it with 
the help of (Q) from the CS fermions to the electrons (CS 
retransformation) . Denoting the expectation value with 
respect to the electronic ground state by (-) e we obtain 



(9) 



m b (J(q,t)) e e z x —(: p(q" + q)p(-q') 

SVO q 



J(q, t) is the electron current operator J(q, t) = 
Efe(fcM)a+ fe> _. /2 (t)a e) g_ g72 (i)+(l/m 6 ) £ f , A(q')p(q- 

q',t). Here A(q') is the external vector potential of the 
v = 1/2 system. Now we use (||) to calculate the first 

term of the right hand side of (||) . Since (J(f, t))Y = = 
for all electron densities of the system one gets zero for 
this term. The cumulant expectation value of the second 
term of ([)]) is of order O(q )(p(q)) . Thus, the third term 
of (Q) is given (q — > 0) by the non-cumulant part 



ipo- 



mb 



-x(f,f)> = -poV 2 (q)(p(q,t)) (10) 



(p(£*)> = 



(v (q) + O(q°)) 



V e *(q,t) 



(11) 



Thus the static density-density response function van- 
ishes as i>o(<?) -1 , i.e linearly for q — > 0. 

The above results are rigorous consequences of the 
equation of motion in the limit q — > 0. We consider 
this result as a consistency condition for an approxima- 
tive calculation of the correlation functions: An approx- 
imative density-density correlation function and the cor- 
responding density-momentum correlation function have 
to satisfy the above analysis. Especially we have seen 
that the 1 /q 2 singularity of the non cumulant part of the 
commutator in (|7j) , due to the commutation with the 
v 2 -tcrm of H (second term in (g)), is cancelled by the 
1/q singularity of the non cumulant part, due to the 
commutation with the ui-term of H (third term in (|j)). 
Thus there is no l/q 2 -singularity in the denominator of 
( |TT| ) . It is clear from the structure of the operators and 
from the derivation above that the first term gets this 
1 /q 2 singularity by averaging with respect to every state 
not only the CS ground state. This is no longer true for 
the second term. In this term one gets the 1/q 2 singular- 
ity by averaging over the CS ground state. This ground 
state is reached by the dynamics of the CS system. Thus 
the consistency relation is written for q — > as 



q'q 
q 



w 2(<?')( : p(q\t)p(q- q'>t) 
2? m b \q ) \ q 

-E— < : (^x*(q + q',t))p(-q',t) :) 

f^ m b \q J V q 



0{q) ■ 



(12) 



Here, the third term cancels the cumulant part of the 
second term to order O(q). 

We add a remark concerning the Hamiltonian (Q). The 
Hamiltonian (0) is truncated by keeping only terms up 
to second order in the quadratic density fluctuations. We 
can show by the same methods as above that an analysis 
of the full problem with the Hamiltonian (S) leads to the 
same results. 

In the next section we formulate a Singwi-Sjolander 
theory of the half-filled Landau level, where we make 
an approximation of the momentum-density and density- 
density correlation function which respects the consis- 



tency condition (12) derived above 
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III. CS RESPONSE FUNCTION WHICH 
INCLUDES DYNAMICAL CORRELATIONS 

In this section, we will calculate response functions of 
the v = 1/2 CS system which include correlations be- 
yond the RPA. As mentioned earlier, one transforms the 
Hamiltonian of electrons in a magnetic field to a CS- 
Hamiltonian (Ji|) at zero magnetic field. So, one can cal- 
culate response functions of the v = 1/2 system with 
approximative methods which were developed for the 
Coulomb system earlier. In this section, we will ap- 



ply the theory of STLS (12| to the CS system in the 
spirit of CC |11[ . The response function matrix \ relates 
the density p(k,uj) and transverse momentum response 
7tt(/c,w) = k/k x jf(k,uj) to an external perturbation, a 
scalar potential V cxt and a transverse vector potential 
via 



p(k,u) 
■K T (k,ui) 



V cxt (k,Lu) 
Af\k,uj) 



(13) 



Following the original derivation of STLS jl2j we start 
from the equation of motion for the one-body Wigner 
distribution function 

k 

which determines the density p(f,t) — ^2ff^{r,$t) 
and the momentum 7?(r, t) = YlpPJ (fiPtt)- The 
Heisenberg equation of motion for f^(f,p;t) is 



^) { f,p-t) = P —f {X) imt) 
at nib 



(15) 



+ / d 2 r 



p' 



(P-P')i 
nib 



Pj_ 



+ V P Af ll) V r ,iV™\?,t) + ^V p ,ifWv r ,iAf*(r,t), 



where V p ,j = d/dpi, and 



k.k' 



x a + -. , (t)a~, . t, ,Jt) :) 

p'-k'/2 y ' p'+k'/2^ I I 



(16) 



is the two-body distribution function. 
Now we have to decouple the correlation function 
fV)(f,p-r',p'-t) in (H). For the Coulomb theory STLS 
uses the following decoupling 

fW(r,p;ripi;t) » f^(r,P;t)f (1) (r',p';t)g(f- r') . 

(17) 

If one uses the decoupling function g(r — f') = 1 in ( |l7j ) 
it is easily seen that one gets with the help of this decou- 
pling the density-density response function in RPA from 



the equation of motion ( |l5| ) . Such a decoupling does not 
respect the small distance correlations. STLS take into 
account these correlations by taking for g(f—r') the equi- 
librium, static pair correlation function. This could also 
be motivated as follows. One can determine g(r — f') by 
summing jl7|) over p, and p' such that both sides of the 
approximation (fiT 



coincide for V oxt — A^ = 



W))W')) 



(18) 



The ansatz (|17|), ( |l8j ) is a decoupling, specific for the 
Coulomb theory because the interaction of the Coulomb 
Hamiltonian consists only of a density-density vertex. In 
the CS theory the Hamiltonian ([|) has in addition to 
the density-density vertex a density-momentum vertex. 
The effect of this is given by the first term in the square 
brackets on the right hand side of equation ([l5|). It is 
not clear whether the decoupling (|l7j), (|l8| ) is a good 
approximation for this term. In their approximation CC 
glj used the decoupling |L7]), |[|) for the CS theory. 

In the following, we want to check if this decoupling 
respects the consistency condition (|l2|). For this, we have 
to determine the density-density correlation function and 
the density-momentum correlation function from the de- 
coupled two body Wigner distribution function. By a 
summation over p, p' we get from (17) 



{:p(f',t)p(r,t) :} 
t)p{r,t) :) 



(p(r',t))(p(r,t))g(r'-r),(19) 
(Tf(r',t))(p(r,t))g(r' -r). (20) 



By inserting these decouplings into (|T^) we get for the 
first summand poV2(q)q(p(q,t)) + 0(q). The third sum- 
mand is zero. The second summand of (|12) is given by 
a similar discussion as in the last section (see especially 
the discussion below equation (^|)) 

. \ - vi(q + k) ( jq + k) \fQfu \ 



m b y \q + k\ 

PoV2(q)q(p(q,t)) + 0(q). 



(21) 



where S(k) — 1 is the Fourier transformion of g(r) — 1. 

With the help of J2ki S ( k ) ~ *) = ~ l we § et for @ 
— (l/2)p V2{q)q(p(q 7 t)) +0(q). Thus we obtain that the 
decoupling of CC does not fulfill the consistency con- 
dition (|I2"|). That is the reason why CC obtained in 
the scheme of their decoupling a zero-frequency density- 
density response function which vanishes as the square 
of the wave vector in the long wavelength limit. 

To get a better insight into the violation of the con- 
sistency condition (|l2|), we calculate the Fourier trans- 
formation of the first two terms in (12) multiplied by q. 
This is given by 



d 2 r'V rii {(:(p(f", t) - po) p{r, t):) V rii v 2 (r ' - r)} 



(22) 



- d 2 r'V r A (:(Vw x fi(f t)) p{r, *):) 



^ r ,iV2(r' - r) 
4irp 
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Because of the appearance of the 7? in the second term in 
( f22| ) we proceed similar as in the discussion below ([)]) by 
retransforming to the electronic system. Then, the term 
of quadratic density fluctuations cancels the first term 
in (p2|). This is no longer true by using the decouplings 
(|l9| ), (p0| ) in (|22|). The operator V r < in the second term 
in (|22|) acts in this decoupling on (^(r'))(p(f)) as well 



as on g(f 



This results in two summands. When 



we take the quadratic density fluctuation part after a CS 
retransformation of these two terms we get that the first 
term is cancelled by the first term in (p2|) (with the help of 



X(f>(r' — r) = 25(r'—r)). The term with V r ' acting on 
' — r) is the reason for the violation of the consistency 



g(r 

condition. In other word doing the CS-retransformation 
first and then the decoupling of the correlation function, 
or vice versa, leads to different results. We require that 
these two actions commute and that fixes the density- 
momentum decoupling in the case of a given decoupling 
of the density-density correlation function. For the given 
density-density decoupling of STLS (|l^) we get by this 
requirement the decoupling 



(:V r / x ir(f" ,t)p(f,t):)i 



Thus we get that g(f' — r) should not be differentiated. 
In the decoupling of the equation of motion (|l^) , we get 
for the second term in the square brackets after a Fourier 
transformation with respect to r: 

J d 2 rd 2 r' e&\y r M?-? , Wp,if ( -* ) (f,Pii f ',P , \t) 

,/,w/v>«nv r ,,^- r -<,] 



(24) 



where v = vq + v 2 ■ The first term in the square brackets 
should be decoupled by 



d 2 r d 2 r' Y (V r 4 & ) (r - r ') 

v 

^ 1 



(25) 



47T j O 



= J d 2 rd 2 r'J2 e W (V r ,iV 2 (f-f')) 
p' 

• (V r x p + V r / xp' + iqxpj W p ^f {2) (f,p;f',p';t) 
« / d 2 rd 2 r'Y e l ^{V rl v 2 (r-r'))-^—g{r-r') 

J 47TP0 

■{V r xp + V rl xp' + iqxp^pJ {1 \r,p-t)f^\r'J'-t), 

where g(f ' — r') is not differentiated. 

With the help of the Heisenberg equation of motion 
( |l5| ) and the approximations (^3) and (^5|) it is possi- 
ble to calculate the response matrix (O) (similar to the 



calculations of CC pdj). By doing this we get 



<g(r' -r)((V r , xn(r>,t)))(p(r,t)). 

(23) 



x = x° [i - u x °] 



where 



(26) 



(27) 



X°p P is the ideal gas density-density response and Xt * s 
the corresponding transversal momentum-momentum re- 
sponse. These response functions are known analytically 
. The matrix of the effective potentials is 



U : 



wo(k) + w 2 {k) iwi(k) 



-iw\{k) 







(28) 



where w a (k) = (1 — G a (k))v a (k) and the local field fac- 
tors G a (k) are given by 

c , 11) 4 E , 1 - sw ^|gf . m 

PQp k a "\k - p\ b <* 

Here ao = 6o = lj Oi = 0, &i = 2 and a 2 = 0, b 2 = 2. 

In the following, we denote the quantities which are 
calculated by the means of CC with an upper index CC, 
for comparison. The structure factor S(k) can be calcu- 
lated with the help of the fluctuation-dissipation theorem 



1 f°° 

S(k) = / dujlm[x PP (k,uj)}, 

Pqk Jo 



(30) 



where x^p is density-density part of the response matrix 
in 0) 



Xin® 



X°pp(k, ijj) 



Xpp(kl U) 1 - X°p P (wo(k) + w 2 (k) + wi (fc) 2 Xr) ' 



(31) 



By using their decou plin g, CC calculated a similar re- 
sponse as in d|) and ((29) with afi c = bfi c = 1, af c = 1, 



6f c = 1 and = 0, 62 = 2. Using the rotational in- 
variance of S(k) it is easy to show that Go(k) = Gq C (A:) 
is linear in k. Further Gi(fc) = G^ik) = G 2 (fc) is 
quadratic in k and Gf c (0) = 1/2. Thus we obtain for 
our decoupling a zero-frequency density-density response 
function which vanishes linearly in the long wavelength 
limit. For the CC decoupling it vanishes as the square 
of the wave vector in the long wavelength limit. Since 
the only difference in the G-terms of our decoupling and 
the decoupling of CC is given by Gi, we will compare 
Gi for these two decouplings. One gets from the integral 
d|) that Gi(fc) approaches to Gp c (fc) for k > k F . Both 
functions are equal for k — > 00. As mentioned above this 
is not the case for k — > 0. 

In the following, we will calculate a solution of the 
(integral) equations (|29|-pT|) by a numerical iteration 
method. As a starting point we use G a = for all a, cor- 
responding to the RPA. For doing this one has to choose 



the dimensionless coupling strength r s = l/(a 0y /np Q ), 
where ciq is the Bohr radius ao = e 2 m&. Since the results 
show little variation with this coupling strength (in the 
region 1 < r s < 10) we will choose for dehniteness r s = 6 
in the following figures. This special choice of r s could 
be motivated through calculations of the effective mass of 
the u =1/2 system |n| . In the following we present and 
discuss our results for various quantities derived from the 
density-density response function xpp- In Fig. [I] we show 
(irpa)S(k, u>), where 



S(k,u) = -— lm[ XpP {k,Lj)} (32) 



is the dynamical structure factor. With the help of 
X° pp = Pak 2 / (m b uj 2 ) and \t = 2-k p\k 2 / {m\u} 2 ) which is 
valid for k/kp <C 1 , kkp/mt, <C oj, we get from ( |3~ij ) that 
S(k,u) has a pole at the cyclotron frequency u = k F /mb 
describing the inter Landau level excitation, which is un- 
affected at k = by correlations and is in agreement with 
Kohn's theorem. For finite values of k we get from Fig. || 
that the cyclotron mode (the pole of S(k,u>)) of CC and 
our decoupling is given by a smaller frequency than in the 
IIP A. From Fig. | we obtain that S(k,w) « S cc (k,uj) for 
lj kkp/nib. This can be understood by the asymptotic 
form of and Xt m this range and the similarity of 
Go(fc) and G2(k) for CC and our decoupling for k <C kp. 
It is also seen from Fig. | that S(k,u) pb S RFA {k,iv) for 
small values of io2m/k 2 F . This can be also seen in Fig. 

\ where we show the function Xpp(fc>0). One sees from 
this figure that x pp (fc,0) w Xp P PA (fc,0) for k < kp and 
Xpp (k,0) w Xpp(k,0) for fc > kp. This is understandable 
by the asymptotic forms of the G(k)s (see the discussion 
above). One also obtains from this figure that x cc (fc,0) 
has the asymptotic 0(k 2 ) for fc — ► 0. 
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FIG. 1: The dynamic structure factor S(k,u>) times Trpo for 
k = 0.6/cf and rs = 6, as a function of uj2m/k 2 F for RPA, 
CC and our decoupling (solid line). The 5- function peak cor- 
responding to the inter Landau level mode has been artifi- 
cially broadened for clarity and contains most of the spectral 
strength. 



TABLE I. The Coulomb interaction energy (e Coul (r s )) per 
particle 



r s 


This Work 


CC RPA 





-1.00e 2 /(aor s ) 


-1.00e 2 /(a r a ) -1.19 t 2 /(a r s ) 


6 


-1.04e 2 /(oor s ) 


-1.04 6 2 /(aor s ) -1.76 e 2 /(a r a ) 



By using the compressibility sum rule @ ^ 

X PP (k,0) \ 



K = =■ lim 

p 2 k^O 



l + 2ir% XpP (k,0) 



(33) 



we can calculate the compressibility K of the v = 1/2 sys- 
tem. The denominator is due to the fact that by changing 
the area of the system the positive background has to be 
changed also to conserve neutrality. With the help of 
Xpp(M) = -m 6 /(27r), Xr(M) - ~Po/m b and © this 
leads to 



K = 



1 



2tt P 2 



2 16 
3 



2tt Po 



oCoul 



(r.)) 



(34) 



We call (e Coul (^)) = l/2^2 n v [k)(S{k) - 1) the 
Coulomb interaction energy per particle, following p9[ . 
For a comparison of the compressibility with the other 
theories we mention that K KPA = (3/4)(m;,/27rpQ) for 
the RPA and K GC = for CC. Now we calculate 
(e Coul (r s )) for RPA, CC and our decoupling. The 
Coulomb interaction energies for r s — ► and r s = 6 
are shown in Table I. We see from this table that the 
Coulomb interaction energies of our theory and CC are 
equal. By calculating the full Coulomb energy per 
particle via coupling constant integration u Coul (r s ) = 
l/(a r 2 ) J ra dr' s (a r' s )(e Coul {r' s )) we get for CC and our 
theory u Coul (r s ) » -1.02 e 2 /{a r s ). The Coulomb en- 
ergy of electrons in the lowest Landau level was cal- 
culated earlier by Morf and d'Ambrumenil jl3| and 
by Girlich E4] by numerical diagonalization methods. 
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FIG. 2: The cyclotron peak of RPA, CC and our decoupling 
(solid line) for r a = 6. It is determined through the singularity 
of Im[x P p]. 
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FIG. 3: Static density-density response function X PP {k, 0) as a 
function of k/kp for the RPA, CC and our (solid line) decou- 
pling scheme for r 3 — 6. The discontinuity in the derivative 
at k = 2kp is due to the Fermi surface. 



Within this method the Coulomb energy per electron is 
given by u^°^(r s ) « — 0. 88e 2 / (a r s ). This energy has to 
be compared with the r s — > Coulomb energy of RPA, 
CC and our method. We obtain that the r s — > energy of 
the STLS-type methods is in a better agreement with the 
numerical results than the Coulomb energy of the RPA. 
(The formula for the coupling constant integration leads 
for small r s to the equality of the Coulomb interaction 
energy with the Coulomb energy.) Furthermore, we see 
from table I that in the STLS-type theories the Coulomb 
energy of the lowest Landau level is a very good approx- 
imation to the total Coulomb energy including higher 
Landau levels. This is not the case for the RPA. To get 
the reason for the equality of the Coloumb energies of 
our method and CC we show in Fig. |J for r s — 6 that 
S(k) = S cc (k) for almost all k/k F . We verify that all 
curves in the figure obey the leading (k/kp) 2 /2 behaviour 
for small k/kF, required by the Kohn theorem. Further- 
more, we find for our method that the w-momenta of the 
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FIG. 4: Static structure factor S(k) as a function of k/kp in 
RPA, CC and our (solid line) scheme for r s — 6. 



dynamical structure factor S{k,Lj) for small k (exclud- 
ing the cyclotron contribution) coincide with that of the 
RPA @, [TJ because for small k our effective potentials 
coincide with that of the RPA. 

At last we calculate the Landau energy. As in the 
case of the calculation of the Landau energy in RPA [Q , 
one obtains with the help of a coupling constant integra- 
tion for the Landau energy in our decoupling scheme an 
infrared diverging energy. This is caused by the simplifi- 
cation in considering only quadratic density fluctuations 
in the derivation of the Hamiltonian (Q) from the exact 
CS Hamiltonian (||) (compare the remark at the end of 
section II). The Hamiltonian has two defects: (i) The 
restriction to the quadratic density fluctuations leads to 
an ultraviolet divergence, (ii) The order of the operators 
in the product of the two p is changed to normal order 
which eliminates the ultraviolet divergence but results in 
an infrared divergence [§). A solution of this problem 
for the energy in the RPA was given in [^| by consider- 
ing for the energy calculation all maximal divergent dia- 
grams (RPA) together with all first order diagrams of the 
Hamiltonian (Q). The first order energy (per particle) of 
the CS Hamiltonian (|) is @ uf mag = 4.00/(m h (a r s ) 2 ). 
In the following, we make a similar calculation for our 
decoupling scheme. Thus, we have to calculate the Lan- 
dau energy per particle u mag = uf mag + ud mag of the 
Hamiltonian (|^) but to extent the first order terms to 
those of the Hamiltonian (H) with the full three parti- 
cle interaction. To this end we have to multiply Ui(fc) 
and v 2 (k) in ( p8|) by a parameter A. Then the solu- 
tion for x °f the (integral) equations ( p9[j3l| ) depends 
on A. From the Hamiltonian (|) we get that for deter- 
mining ud mag one has to calculate two terms (one con- 
tains the density-density response function x PP (/c,w;A), 
the other contains the density-momentum reponse func- 
tion Xpir(fc, A)). After some algebra one gets for ud mag 
through coupling constant integration 

ud mag = 



1 OO 

k 



-v 2 {k) Im[x P p(fc,w; A) 



X° PP (k,u)] 



v 1 (k)w 1 (k) Im[xT(k,u)x PP (k,uj; A)] 



(35) 



With the help of a numerical integration of Q35| ) we get 
ud mag = -2.16/(m b (a r s ) 2 ). The total Landau energy 
per particle u mag is then given by 



uf raag + ud 1 



1.84 



1 



m b (a r s ) 2 



(36) 



By doing a similar calculation for the total Landau en- 
ergy per particle within the RPA || we get u^p^ = 
1. 60 /(mb{aor s ) 2 ). The Landau energy can be compared 
to the exact Landau energy oj c /2 of the v = 1/2 sys- 
tem. This is given by u™ ag = 2.00/{m b (a r s ) 2 ). Thus 
we obtain that, in comparison with the RPA, the Lan- 
dau energy and the Coulomb energy of the STLS-type 



8 



theories are in better agreement with the exact Landau 
energy and the Coulomb energy (from numerical diago- 
nalization) . 



IV. CONCLUSION 

In this paper, we consider an approximation of the re- 
sponse function of the v =1/2 system which goes beyond 
the RPA. The method we used is the STLS-theory, first 
established for the Coulomb system fl2|| . Recently this 
theory was applied to the CS system by CC [jll). In 
their theory, CC obtain a density-density response func- 
tion which vanishes as the square of the wave vector in 
the long wavelength limit. We show in this paper that 
the zero- frequency density-density response function van- 
ishes linearly in the long wavelength limit independent of 
any approximation. We obtain this result by considering 
the Heisenberg equation of motion for the longitudinal 
momentum operator. From this equation of motion, we 
derive a consistency condition for a decoupling of the 
equal time density-density and density-momentum corre- 
lation functions. We show that this consistency condition 
is not fulfilled in the theory of CC and that is the reason 
for the quadratic behaviour of the zero-frequency long 
wavelength limit of the density-density response. Based 
on the functional form of the Heisenberg equation of mo- 
tion of the Wigner distribution function (with external 
potentials), we suggest a decoupling of the correlation 
functions in this equation which respects the consistency 
condition. 

We solve the decoupled Heisenberg equation of motion 
by numerical iterations and get the response functions of 
the theory. In contrast to the theory of CC, we obtain a 
density-density response function which vanishes linearly 
in the long wavelength limit for zero-frequency. We get 
agreement for the density-density response function with 
the theory of CC for momenta k > kp. For k < kp we 



get agreement with the density-density response function 
of CC for frequencies u> > kkp/rrn, and to the RPA for 
lu < kkp/nib- Further, we calculate the compressibil- 
ity of the theory by using the compressibility sum rule. 
We obtain a Coulomb correction to the compressibility 
not contained in the RPA. With the help of the response 
functions, we calculate the static structure factor, the ex- 
citation spectrum and the Landau as well as the Coulomb 
energies. As in the theory of CC, we get density excita- 
tions which are lower in their frequencies as a function 
of the wave vector than the excitations calculated with 
the help of the RPA. The obtained excitation spectrum is 
almost identical to the spectrum of CC. The same holds 
for the static structure factor. This is the reason for the 
agreement of the lowest Landau level Coulomb energy as 
well as the full Coulomb energy with the energies calcu- 
lated by CC. The lowest Landau level Coulomb energy 
fits better to the Coulomb energy calculated by numerical 
methods (jJHH (w 114% of the numerical lowest Lan- 
dau level Coulomb energy) than the RPA. We remark 
that the relative part of the Coulomb energy in the low- 
est Landau level, i.e. linear in e 2 , is much enhanced as 
compared to the RPA. Finally, we calculate the Landau 
energy of the system. We obtain a much better approxi- 
mation of the exact Landau energy ( ss 92% of the exact 
Landau energy of the v = 1/2 system) than the RPA. 

In summary, a consistent decoupling of the Wigner 
function in the Heisenberg equation of motion leads to 
results which are in a better agreement with known nu- 
merical and exact results, respectively, than the RPA. 
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